The chemistry of La on the Si(OOl) surface from first principles 
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This paper reports state-of-the-art electronic structure calculations of La adsorption on the Si(001) 
surface. We predict La chains in the low coverage limit, which condense in a stable phase at a 
coverage of | monolayer. At | monolayer we predict a chemically rather inert, stable phase. La 
changes its oxidation state from La 3+ at lower coverages to La 2 + at coverages beyond 1 monolayer. 
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I. INTRODUCTION 

Device scaling has been the engine driving the micro- 
electronics revolution as predicted by Moore's lawi By 
reducing the size of transistors, processors become faster 
and more power efficient at an exponential rate. Cur- 
rently the main challenge in device scaling is the integra- 
tion of high-k oxides as gate oxides into silicon technol- 
ogy 

The gate oxide is an integral part of a metal-oxide- 
semiconductor field-effect transistor (MOSFET). It is the 
dielectric of a capacitor, which is used to attract charge 
carriers into the channel between source and drain, and 
thus switches the transistor between its conducting and 
its non-conducting state. With a thickness of approxi- 
mately 1-2 mn^ the gate oxide is the smallest structure 
of a transistor. Further scaling would result in an unac- 
ceptably high quantum mechanical leakage current and 
thus a large power consumption. 

In current transistors, the gate oxide is made from Si02 
and SiO^Nj,. Future transistor generations will have to 
employ oxides with a higher dielectric constant (high- 
k). This allows greater physical thicknesses and thus 
reduces the quantum mechanical leakage currents. The 
main contenders for the replacement of Si02 in future 
transistors are, from today's point of view, oxides con- 
taining alkaline earth metals like Sr or Ba, third-row el- 
ements like Y or La, forth-row elements like Ti, Zr and 
Hf, or mixtures thereof. Prominent examples are per- 
ovskite structures around SrTiOg^ and LaAlOg 4 * 5 ^, flu- 
orite structures like Zr02 and HfO^i and also Y2O3 and 
La20g^ or pyrochlore structures like La2Hf207^£ and 
La2Zr2C>7iii*i£ Recently, also promising results on P^O,-? 
have been published^ 

While the first high-k-oxides will be grown with an in- 
tcrfacial Si02 layer, a further reduction in scale requires 
high-k-oxides with a direct interface to silicon. The re- 
quirement to limit interface states, and the often crys- 
talline nature of the oxides demand an epitaxial growth of 
the oxides on silicon. Considering layer-by-layer growth 
by molecular beam epitaxy (MBE), the first growth step 



for high-k oxides is the deposition of the metal on sili- 
con. Therefore we have investigated deposition of metals 
out of the three most relevant classes for high-k oxides 
on Si(OOl). These are the divalent alkaline-earth met- 
als and the three- and the four-valent transition metals. 
The results on adsorption of Zr and Sr have been pub- 
lished previouslyii 4 * 1 ^ The present paper completes the 
study with a description of La-adsorption on Si(OOl) as 
example of a trivalent metal. 

Our previous work has shown that Zr tends to form 
silicides readily. 14 Silicide grains have been observed af- 
ter Zr sputtering on Si(001)r& unless silicide formation 
is suppressed by early oxidation which, however, leads to 
interfacial SiC>2. The Sr silicides are less stable in con- 
tact with silicon and due to their sizable mismatch in 
lattice constant, nucleation does not proceed easily. The 
alkaline-earth metals Sr and Ba have been used in the 
first demonstration of an atomically defined interface be- 
tween a high-k oxide, namely SrTiC"3 and silicon^ By 
following through the detailed steps of the formation of 
this interface, starting at the low-coverage structures of 
metal adsorption, we were able to provide a new pic- 
ture for the phase diagram of Sr on Si(001)J^ The phase 
diagram has been important to link the theoretical inter- 
face structure of SrTiC>3 on Si(OOl) to the experimental 
growth processi 1 ^* 1 ^* 1 ^ From the interface structure and 
its chemistry we could show in turn that the band-offset, 
a critical parameter for a transistor, can be engineered to 
match technological requirements by carefully controlling 
the oxidation of the interface H 

Since many of the characteristics of Sr adsorption carry 
over to La-adsorption let us briefly summarize the main 
results. 15 Sr donates its electrons to the empty dan- 
gling bonds of the Si-surface. The Si-dimers receive 
electron pairs one-by-one, and unbuckle as they become 
charged. When all Si dangling bonds are filled, i.e. be- 
yond i monolayer (ML), additional electrons enter the 
anti-bonding states of the Si-dimers at the surface, and 
thus break up the Si-dimer-row reconstruction. 

At low coverage, Sr forms chains running at an angle 
of 63° to the Si-dimer rows. As the coverage increases, 



2 



the chains condense first into structures at g ML and at 
j ML, which are determined by the buckling of the Si- 
dimers and their electrostatic interaction with the pos- 
itive Sr ions. At \ ML a chemically fairly inert layer 
forms, where all dangling bonds are filled and all ideal 
adsorption sites in the valley between the Si-dimer rows 
are occupied. 

The paper is organized similar to our previous work 
on Sr adsorption. In Sec. ^ we describe the computa- 
tional details of the calculation. In Sec. IIIII and IIVI wc 
review briefly the reconstruction of the Si(OOl) surface 
and we discuss the known bulk La silicides. Sec. ]V\ IV1I 
and IVHI deal with the low coverage limit, where La ad- 
atoms form dimers and chain structures. Beyond the 
canonical coverage of 1/3 ML fSec. IVIIl)) we observe a 
change in the oxidation state of the La ad-atoms from 
3+ to 2+ (Sec. IIX|) . The results are placed into context 
in Sec. Q(] where we propose a phase diagram for La on 
the surface. The computational supercells used for the 
simulation of the low-coverage structures are shown in 
the appendix. 



II. COMPUTATIONAL DETAILS 

The calculations are based on density functional 
theorjii2iS£ using a gradient corrected functional i2I The 
electronic structure problem was solved with the projec- 
tor augmented wave (PAW) method, 22,23 an all-electron 
electronic structure method using a basis set of plane 
waves augmented with partial waves that incorporate the 
correct nodal structure. The frozen core states were im- 
ported from the isolated atom. For the silicon atoms we 
used a set with two projector functions per angular mo- 
mentum for s and p-character and one projector per an- 
gular momentum with d-character. The hydrogen atoms 
saturating the back surface had only one s-type projector 
function. For lanthanum we treated the 5s and 5p core 
shells as valence electrons. We used two projector func- 
tions per magnetic quantum number for the s, p, and / 
angular momentum channels and one for the d channel. 
The augmentation charge density has been expanded in 
spherical harmonics up to £ = 2. The kinetic energy cut- 
off for the plane wave part of the wave functions was set 
to 30 Ry and that for the electron density to 60 Ry. 

A slab of five silicon layers was used as silicon sub- 
strate. This thickness was found to be sufficient in pre- 
vious studies on Sr adsorption^ The dangling bonds of 
the unreconstructed back surface of the slab have been 
saturated by hydrogen atoms. The lateral lattice con- 
stant was chosen as the experimental lattice constant 
a = 5.4307 A of silicon, 24 which is 1 % smaller than 
the theoretical lattice constant. Since we always report 
energies of adsorbate structures relative to the energy of 
a slab of the clean silicon surface, the lateral strain due 
to the use of the experimental lattice constant cancels 
out. The slabs repeat every 16 A perpendicular to the 
surface, which results in a vacuum region of 9.5 A for the 



clean silicon surface. 

The Car-Parrinello ab-initio molecular dynamics^ 
scheme with damped motion was used to optimize the 
electronic and atomic structures. All structures were 
fully relaxed without symmetry constraints. The atomic 
positions of the backplane of the slab and the terminating 
hydrogen atoms were frozen. 

Many of the adsorption structures are metallic, which 
requires a sufficiently fine grid in k-space. We used an 
equivalent to twelve by twelve points per (lxl) surface 
unit cell. Previous studies have shown that a mesh of 
eight by eight k-points is sufficient ~ We have chosen a 
higher density here as this allows us to use commensurate 
k-meshes for 3x and 2x surface reconstructions. 

For metallic systems, the orbital occupations were de- 
termined using the Mermin functional, 26 which produces 
a Fermi-distribution for the electrons in its ground state. 
The electron temperature was set to 1000 K. In our case 
this temperature should not be considered as a physi- 
cal temperature but rather as a broadening scheme for 
the states obtained with a discrete set of k-points. The 
Mermin functional adds an entropic term to the total 
energy, which is approximately canceled by taking the 
mean of the total energy U(T) and the Mermin-free en- 
ergy F(T) = U(T) - TS(T) as proposed by GillaniSi 



U(T = 0)n~(F{T) + U(T)). (1) 

The forces are, however, derived from the free- 
energy F(T). The resulting error has been discussed 
previously^ 

In order to express our energies in a comprehensible 
manner, we report all energies relative to a set of ref- 
erence energies. This set is defined by bulk silicon and 
the lowest energy silicide LaSi2. The reference energies 
are listed in Tab. [I] The reference energy .Eg [La] for a 
La atom, corresponding to the coexistence of bulk silicon 
and bulk La, is extracted from the energy E[LaSi2] of 
the disilicide calculated with a 9 x 9 x 3 k-mesh for the 
tetragonal unitccll with a = 4.326 and c = 13.840 and 
the reference energy of bulk silicon Eq [Si] as 

£b[La] = £[LaSi 2 ] - 2E [Si]. (2) 

The bulk calculation for silicon was performed in the 
two atom unit cell with a (10 x 10 x 10) k-mesh and at 
the experimental lattice constant of 5.4307 Am^ 

For the surface calculations, we always subtracted the 
energy of a slab of the clean (4 x 2) silicon surface of the 
same layer thickness and backplane. 

III. THE SILICON SURFACE 

Before discussing the adsorption of La, let us briefly 
summarize the chemistry of the clean (001) surface 
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TABLE I: Reference energies used in this paper without 
frozen core energy. See text for details of the calculation. 



TABLE II: Energies per La atom of bulk silicides relative to 
our reference energies (Tab.^J. 





Energy [H] 




E[La] [eV] 


So [Si] 


-4.0036 


LaSi 2 (141/amdS) 


0.00 


So [La] 


-32.1395 


LaSi (Pnma) 


0.42 


£ ( ( ) lxl) [5 layer-Si-slab] 


-21.1139 


LaSi (Cmcm) 


0.62 



La 3 Si 2 (P4/mbm) 



0.80 




FIG. 1: Top view of the Si(001) surface and the four high 
symmetry positions spanning the surface irreducible (2x1) 
unit cell. The dimer buckling is not shown. 



of silicon. A more detailed account has been given 
previously^ 

On the unreconstructed silicon surface, the atoms form 
a square array. Due to a lack of upper bonding partners, 
each atom has two singly occupied dangling bonds point- 
ing out of the surface. In order to avoid half-occupied 
bands, pairs of silicon atoms dimerize, using up one dan- 
gling bond per atom to form the dimer bond. This 
is called the dimer row reconstruction. Still, one dan- 
gling bond per silicon atom is half occupied, which drives 
the so-called buckled-dimer reconstruction: One atom of 
each dimer lifts up and the other shifts down, resulting in 
a "buckled" dimer. This buckling causes both electrons 
to localize in the upper, sp 3 -like silicon atom of a dimer, 
whereas the other, sp 2 -like silicon atom with the empty 
p-like dangling bond prefers a more planar arrangement. 

Fig. ^ shows a ball stick model of the silicon surface 
and introduces the high-symmetry adsorption sites A to 
D, which we will refer to in the following discussion. 

Filling the empty dangling bond with two electrons 
results in a removal of the buckling as observed in our 
studies on Sr adsorption^ As La has an odd number of 
valence electrons, we also examined the changes of the 
buckling upon filling the initially empty dangling bond 
with a single electron. While the difference in z coordi- 
nate of the two silicon atoms of a buckled dimer is 0.76 A 
and the one of an unbuckled dimer 0.00 A, it is 0.35 A 
after donation of a single electron to a dimer. Thus the 
amplitude of the dimer buckling may be used as a mea- 
sure for the electron count. 



IV. BULK LA SILICIDES 

In the case of Sr, the chemical interaction with sil- 
icon could be understood by investigating the bulk Sr 
suicides^ All these structures could be explained by 
the Zintl-Klemm concept^ The electropositive Sr atoms 
donate their two valence electrons to the silicon atoms. 
Each accepted electron saturates one of silicon's four va- 
lences. A Si~ can thus form three covalent bonds, a Si 2 ~ 
only two, a Si 3- only one and a Si 4- has a closed shell 
and does not form covalent bonds. This principle was 
found to also be valid for the surface reconstructions of 
Sr on Si(001) 

The Zr silicides on the other hand cannot be explained 
by the simple Zintl-Klemm concept The Zr d states 
also contribute to the bonding and thus retain a variable 
number of electrons. 

Similarly, the La silicides cannot be simply explained 
by a quasi-ionic interaction with silicon. We find La in 
formal charge states between two and three (i.e. charge 
according to the Zintl-Klemm concept). Also the atom 
and angular momentum resolved density of states reveals, 
that La d states are partly occupied in these structures. 

Fig.|21shows the La silicides. LaSi2 is the lowest energy 
silicide. The energies per La atom are listed in Table ILTI 



V. AD-ATOMS AND LA DIMERS 

Our search for the adsorption structures of La have 
been guided by the electron count rules that emerged 
from our investigation of Sr— adsorption on the same 
surface. The studies of Sr provided a consistent picture: 
The electrons from Sr are fully transfered into the Si- 
dimer dangling bonds of the Si substrate. The ordering 
of Sr atoms on the surface is determined by the electro- 
static attraction between the Sr-cations and negatively 
charged Si-ions at the surface. The negative Si-ions are 
the raised atoms of buckled Si dimers and the atoms of 
filled, and thus unbuckled, dimers. This picture holds up 
to coverages where all Si-dimers are filled at 1/2 ML. Due 
to the different electron count of La as compared to Sr, 
we expect that the silicon dimers are filled already at a 
coverage of 1/3 ML and secondly we anticipate deviations 
from the above scheme. 

Even though we predict La-chains to be the most stable 
structures in the low coverage limit, we first investigate 







FIG. 2: (Color online) The bulk silicide structures, top left: 
LaSi2 141/amdSSm (lowest energy structure); top right: LaSi 
Pnmaj— bottom left: LaSi Cmcmj— bottom right: La3Si2 



isolated ad-atoms^ and La-dimers in order to provide an 
understanding of the constituents of the more extended 
structures. Chain structures are more stable by 0.34- 
0.40 eV per La atom as compared to isolated ad-atoms. 

Similar to Sr?i£ we find the most stable position of 
an isolated La atom at position A, in the center of four 
Si-dimers (compare Fig. QJ. The position D, B and C 
have energies 0.23 eV, 0.51 eV and 1.70 eV higher than 
position A. A (4 x 4) supercell has been used for these 
calculations. 

The diffusion barrier along the valley is equal to the 
energy difference between sites A and B, namely 0.52 eV, 
the one across the row is 1.31 eV and is estimated by the 
midpoint between the sites A and D. 

The formation of La dimers lowers the energy per ad- 
atom by 0.10-0.18 eV compared to isolated ad-atoms. 
Due to the topology of the Si(001) surface, three differ- 
ent types of La-dimers can be formed: (1) orthogonal to 
the Si dimer rows, (2) parallel to the Si dimer rows and 
(3) diagonal to the Si dimer rows. All three structures 
are shown in Fig. EJ We find that the parallel La-dimer is 
lowest in energy, followed by the orthogonal and diagonal 
La-dimers. All La-dimers lie within a small energy win- 
dow of 0.08 eV. Note, that we only investigated singlet 
states. 

A pair of La atoms has six valence electrons and from 
the lessons learned from Sr adsorption, one would assume 
that three Si-dimers in the vicinity of the La-dimer get 
unbuckled. This is, however, not the case. Only two 
Si-dimers become fully unbuckled. The remaining two 
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FIG. 3: Schematic representation of isolated La-dimers on the 
Si(001) surface. The filled circles represent the La ad-atoms, 
the rectangles represent a filled and therefore unbuckled Si- 
dimers. The triangles represent buckled Si-dimers. The flat 
side of a buckled dimer indicates the upper Si atom with a 
filled dangling bond, whereas the pointed side indicates the 
lower Si atom with the empty dangling bond. Only the Si- 
dimers which are clearly unbuckled have been drawn as rect- 
angles. The partially unbuckled ones are represented by tri- 
angles (see discussion in section 0. The supercells used for 
the total energy calculations are shown in Fig. 1141 



electrons from the La-dimer enter into states that are 
derived from the upper dangling-bond band and which 
have an admixture of the La-c? and / states. 

Usually one can easily distinguish between buckled and 
unbuckled dimers. In the vicinity of La-dimers oriented 
diagonally or orthogonally, however, we also observe Si- 
dimers with an intermediate buckling amplitude. Thus, 
in these cases, the oxidation state of the La atom, namely 
2+ versus 3+ cannot be attributed in a unique manner. 

For the orthogonal and diagonal La-dimers we observe 
a tendency for the La atoms to reduce their distance com- 
pared to staying centered on A sites by 1 - 4%. For 
the parallel La-dimer this effect is opposite and much 
larger. The distance between the two La atoms is 4.11 A, 
compared to 3.84 A between two A sites which amounts 
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to an expansion of 7 %. This ad-atom repulsion within 
one valley has already been observed in case of Sr— and 
explains the formation of isolated chains instead of con- 
densed chains or clusters at low coverages. Nevertheless, 
we find the parallel dimer to be the most stable. 



VI. CHAIN STRUCTURES AT LOW 
COVERAGES 

As we combined pairs of the most favorable site for 
isolated La ad-atoms into dimers, we now search for ways 
to stack the three types of La-dimers together into more 
extended structures. 

We systematically approached linear chain structures. 
Each of the three La-dimer types - parallel to the Si- 
dimer row, orthogonal or diagonal - has been stacked 
such that it shares at least one Si-dimer, so that this Si- 
dimer is next to two Lanthanum atoms from different La- 
dimers. The energetic ordering has been deduced on the 
basis of binding energy per La atom (compare Tab. IHIfl . 
Note that the binding energy per La atom for a given 
chain structure is slightly coverage dependent. In case 
of the double stepped chains (compare Fig. 0^, for the 
structural principle) the adsorption energy varies within 
0.06eV at coverages between 1/10 and 1/5 ML. In order 
to ensure comparability all numbers listed here refer to 
a coverage of 1/6 ML. 

We start our investigation with parallel La-dimers 
shown in Fig. |3Jd, which is the most stable dimer struc- 
ture. The most favorable chain in this class is stacked 
perpendicular to the Si-dimer rows as shown in Fig. 
Its energy lies 0.06 eV per La atom above the lowest en- 
ergy chain structure. 

The most favorable chain made from orthogonal La- 
dimers is shown in Fig. Qb>. It can also been interpreted 
as a variant of a chain of diagonal La-dimers (compare 
Fig-CDc). This is the most favorable chain structure of La 
atoms on Si(001). Its chains run at an angle of about 76° 
to the Si-dimer row. It should be noted that it is equally 
possible to arrange the La-dimers in a zig-zag manner 
as shown in 0Jx The zig-zag chain has not been explic- 
itly calculated. The coexistence of straight and zig-zag 
chains has been found for Sr on Si(001), where the two 
modifications have been shown to be almost degenerate 
in energy^ 

In all low-energy structures each La atom is thus sur- 
rounded by four silicon atoms having filled dangling 
bonds. Three of them are partners of filled Si-dimers 
while one is a buckled Si-dimer with the negative Si atom 
pointing towards the La ad-atom. On the basis of count- 
ing unbuckled Si-dimers, these structures are in a 3+ 
oxidation state. 

The La-chain is the configuration with lowest energy 
in the low coverage limit. The lowest energy chain struc- 
tures are of the order of 0.17 eV per La atom more stable 
than the most favorable isolated La-dimer. At elevated 
temperatures, entropic effects will lead to increasingly 



La-dimer type 


a E[La] [eV] panel 


parallel 


90 


-0.30 
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parallel 


63 


-0.15 
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parallel 


45 


-0.20 
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parallel 


34 


-0.13 
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parallel 





-0.07 
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orthogonal 


90 


-0.26 
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orthogonal/diagonal 


76 


-0.36 


g 


diagonal 


63 


-0.28 


h 



TABLE III: Energies per La atom of the chain structures at 
1/6 ML. The orientation of the chain is described by the angle 
a (degrees) of the chain to the Si-dimer row. The supercells 
used for the total energy calculations are sketched in the cor- 
responding panels of Fig. 1151 
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FIG. 4: Schematic representation of isolated La ad-atom 
chains, a) a single, double-stepped La chain. This is the ener- 
getically most favorable surface reconstruction at low cover- 
ages. A change in chain direction is realized by stacking two 
La-dimers with different orientation (b). Panel (c) shows the 
lowest energy chain structure derived from parallel dimers. 
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the Si-dimer row direction— This restriction has been 
attributed to the requirement that every cation be sur- 
rounded by four Si-atoms with filled dangling bonds, and 
that there is no frustration of the Si-dimer buckling, i.e. 
adjacent Si-dimers are buckled antiparallel. 

For La the situation is more complex. Due to the 
longer periodicity of the La chains compared to those 
of Sr, there are two families of chain packing for La as 
shown in Fig. EJ In the first family the La chains are 
displaced only parallel to the Si-dimer row direction. In 
the second family the chains are in addition displaced 
perpendicular to the Si-dimer row. 

The first family has a preference of (2n + 1) surface 
lattice spacings along the dimer row as in the case of 
Sr adsorption. The spacing in the second family is arbi- 
trary. The reason is that in family one, the buckling of 
every second Si-dimer row is pinned on both sides by two 
neighboring La chains (see Fig.|S^). A Si-dimer is pinned, 
if its buckling is determined by the Coulomb attraction 
of its raised, and thus negatively charged, Si atom to a 
nearby La ion. Since the buckling alternates along the Si- 
dimer row, this pinning can lead to indirect, long-ranged 
interaction between different La-chains. 



FIG. 5: Layer resolved density of states of 1/5 ML. The arrow 
indicates the upper dangling bond bands in the gap of silicon. 
The La panel was magnified by a factor of 2.5. The seemingly 
large gap of the silicon substrate is due to finite size effects and 
also found for the clean silicon surface at this slab thickness. 
For bulk silicon we obtain a "typical" GGA value of 0.65 eV. 
This DOS corresponds to the supercell outlined in Fig. |SJb. 

shorter chain fragments. From the energy-difference of 
the linear chain and the isolated La-dimer, we obtain 
an estimate for the chain termination energy of approxi- 
mately 0.09 eV. It should be noted, that experiments of- 
ten observe shorter chain sequences than predicted from 
thermal equilibrium as high-temperature structures are 
frozen in. 

The electronic structure of the La-chain is analogous 
to that of the Sr single chain— The empty silicon surface 
has an occupied and an un-occupied band formed from 
the dangling bonds of the Si-dimers. La donates electrons 
into the upper dangling bond band. Those dangling bond 
states, which become filled, are shifted down in energy 
due to the change in hybridization on the one side and 
due to the proximity of the positive La-cations on the 
other side. 



VII. CONDENSED CHAINS 

With increasing coverage, the chains become closer 
packed. In the case of Sr, there was a preference for 
a periodicity of (2n + 1) surface lattice spacings along 



In the second family the buckling of every Si-dimer 
row is pinned only at one La-chain as seen in Fig. |SJd, 
while there is no preference of the Si-dimer buckling at 
the other La-chain. Thus for La we find - in contrast to 
Sr— - arbitrary chain spacings. 

The closest packing of La-chains before they merge is 
1/5 ML. We consider two La-chains merged if La atoms 
of different La chains occupy nearest-neighbor A sites 
within one valley. We predict a distinct phase at this 
coverage as seen in Fig. ^] and discussed later. This 
structure, shown in Fig. 03, is derived from chains of 
the second family. An explanation for finding a phase at 
1/5 ML is that the energy at higher coverage increases 
due to the electrostatic interaction of the La atoms within 
one valley. For the first family, the highest possible cov- 
erage before La-chains merge is 1/6 ML (Fig. [B^,). 

Note that the chains can change their direction with- 
out appreciable energy cost as shown in Fig. 0Jd. Exper- 
imentally measured diffraction patterns would reflect a 
configurational average. 

The layer resolved density of states is shown in Fig. 
We see that the Fermi-level lies in a band gap of the 
surface. Above the Fermi-level and still in the band-gap 
of bulk Si, surface bands are formed, which originate from 
the remaining empty dangling bonds of the buckled Si- 
dimers. As in the case of Sr, these states form flat bands 
in the band-gap of silicon, which approximately remain 
at their energetic position as the La coverage is increased. 
Its density of states, however, scales with the number of 
empty dangling bonds. 
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FIG. 6: Schematic representation of the two packing types 
of double stepped La chains at their maximum condensed 
coverage. The reconstruction in panel a) consists of paral- 
lel La chains. The chains in panel b) are also displaced by 
one valley orthogonal to the Si-dimer rows. The ovals indi- 
cate the pinning of the Si-dimer buckling by a La ad-atom. 
The dashed double zig-zag lines shows the positions where, 
in case of structure a), an even number of Si-dimers can be 
inserted in order to arrive at more dilute coverages. In case of 
structure b) an arbitrary number of Si-dimers can be inserted, 
as the buckling of each row is just pinned on one side. The 
calculational supercells are outlined. 



VIII. THE CANONICAL SURFACE AT 1/3 ML 
COVERAGE 



If the spacing of the chains is further reduced, they 
condense at 1 /3 ML to the structure shown in Fig. [7\ 

There are several versions of this structure type. They 
have a repeating sequence of two La-atoms and one va- 
cant A site in each valley in common. The relative dis- 
placement of this sequence from one valley to the next, 
however, may differ. We investigated several structures 
and found the one shown in Fig.0to be the most stable. 

A structure with a sequence of four A sites occupied 
with metal ions separated by two empty A sites, has been 
the most favorable structure at this coverage in the case 
of Sr adsorption. 15 For La, however, this configuration is 




FIG. 7: Schematic representation of the most stable recon- 
struction at the canonical coverage of 1/3 ML. All Si-dimer 
dangling bonds are filled. This structure can be thought of 
as the condensed chain structure in Fig. El b) with a reduced 
chain spacing. The calculational supercell cell is outlined. 



energetically unfavorable. 

At a coverage of 1/3 ML, all silicon dangling bonds 
are filled due to the electrons provided by the La ad- 
atoms. This surface is isoelectronic to the Sr covered sur- 
face at 1/2 MLii 5 - For the Sr-covered silicon surface, the 
increased oxidation resistance of the corresponding 1/2 
ML structure has been observed experimentally. 34 Simi- 
larly we suggest that the surface covered with 1 /3 ML of 
La will have an increased oxidation resistance. 

In Fig. [S] we show the layer-resolved density of states of 
the most stable structure at 1/3 ML. In analogy to the 
1/2 ML covered Sr surface, there are no surface states 
deep in the band gap of silicon, because all Si-dimer dan- 
gling bonds are filled and shifted into the valence band 
due to the electrostatic attraction of the electrons to the 
positive La ions. Note, however, that in contrast to the 
canonical surface coverage of Sr on Si(OOl) at a cover- 
age of 1/2 ML, the canonical La surface exhibits vacant 
A-sites. 



IX. TRANSITION FROM LA 3+ TO LA 2+ 
ABOVE 1/3 ML 

Up to the canonical coverage of 1/3 ML, all thermo- 
dynamically stable reconstructions could be explained by 
La being in the 3+ oxidation state. In contrast to the 
isolated La-atoms and La-dimers, the oxidation state can 
clearly be identified from the number of unbuckled Si- 
dimers: Each unbuckled dimer has received two electrons. 
A 3+ oxidation state is also consistent with the density 
of states. 

If we follow the picture that emerged from Sr, we would 
anticipate that increasing the coverage above 1 /3 ML in 
case of La would lead to filling the Si-dimer antibonds, 
which results in a breaking up of the dimer bonds. For 
La the situation is different: the La-d band is located 
at much lower energies as compared to Sr. Therefore 



a) 



total 






La 




f\ 


top Si 






layer 1 




n ^ 


layer 2 


A A JA 




layer 3 




— . 


layer 4 







-4 -3 -2 
Energy [eV] 



-1 



FIG. 8: Layer resolved density of states of 1/3 ML. The La 
panel was magnified by a factor of 2. Compare Fig. [S] for a 
discussion about the Si band gap. This DOS corresponds to 
the supercell outlined in Fig. |7| 
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FIG. 9: The lowest energy structures at 2/3 ML in the 2+ (a) 
and 3+ (b) regime. The calculational supercells are outlined. 




the energy to break the Si-dimer bonds is larger than 
that to add electrons into the La (i-shell. As a result we 
find that La changes its oxidation state from 3+ to 2+. 
Oxidation states of La that are even lower are unfavorable 
due to the Coulomb repulsion of electrons within the La- 
d and / shells. Thus the structures above 1/3 ML can be 
explained in terms of La 2+ ions and are similar to those 
found for Sr~ 

It may be instructive to compare two structures with 
different oxidation states of La. A good example is found 
at a coverage of 2/3 ML: The lowest energy structure 
is a (3 x 1) reconstruction already found for Sr— and 
depicted in Fig. EJi. This is a clear 2+ structure. Since 
every Si-dimer only accepts two electrons, they can just 
accommodate two of the three valence electrons of La. 
The lowest structure with formal La 3+ ions, which can 
clearly be identified as having all Si-dimer bonds broken, 
is shown in Fig. EJa. It has an energy which is 0.36 eV per 
La atom higher than the structure with La 2+ ions. 

At 1/2 ML, we find a structure where all A sites are 
occupied to be most stable. There the La d-states are 
partially occupied. We confirmed that the system is not 
spin polarized. 

The crossover of the energy surfaces of the 2+ and 
the 3+ structures is shown in Fig. ^3 using a set of sur- 
face reconstructions, for which the charge state can be 



T 




I i I i I i I i I 

0.0 0.2 0.4 0.6 0.8 

coverage [ML] 



FIG. 10: Crossover of the total energy surfaces of the 2+ and 
3+ regime beyond the canonical coverage of 1/3 ML. 



determined unambiguously. It can be clearly seen that 
the 2+ structures become significantly more stable above 
1/2 ML. From Fig. ^] it is apparent that the energy 
rises sharply as the La atoms cross over to the 2+ oxida- 
tion state beyond the canonical interface at a coverage of 
1/3 ML. 
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FIG. 11: (Color online) The surface energy^ versus coverage. 
The open diamonds correspond to the thermodynamically ac- 
cessible surface structures while the triangles mark metastable 
structures. Compare Figs. |SJd,[7| and for the structures at 
1/5, 1/3 and 2/3 ML, respectively. At 0.42 ML we predict 
a (2 x 1) reconstruction which originates from the half-ML 
structural template with a La vacancy concentration of 17 % 
(see discussion in the text) . 
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FIG. 12: The phase diagram for La adsorption on Si(001). 



X. PHASE DIAGRAM 

Based on the surface energies composed in Fig. ITT1 
we extracted the zero-Kelvin phase diagram shown in 
Fig. E| The slope of the line-segments of the lower en- 
velope in Fig. E] corresponds to the chemical potential, 
at which the two neighboring phases coexist (for a more 
elaborate discussion, refer to Ref^). The stable phases 
are defined by the coverages where two line segments with 
different slopes meet. The zero for the La chemical po- 
tential has been chosen as the value at which LaSi2 and 
silicon coexist. Consequently, all phases in regions of 
positive chemical potentials are in a regime where the 
formation of bulk silicides is thermodynamically favor- 
able. 

Below a chemical potential of — 0.44eV we expect sin- 
gle chain structures as described in Sec. IVII At 1/5 ML 
we predict a distinct phase since this is the highest pos- 
sible coverage without La ad-atoms at nearest neigh- 



bor A sites (compare Fig. Et>)- At a chemical potential 
of -0.19 eV the stability region of the 1/3 ML coverage 
(Fig. EI starts. 

The transition from the phase at 1/3 ML to the 
2x1 reconstructed surface at 1/2 ML, where all A- 
sites are filled, can be described by a decrease of La- 
vacancies (compare Fig. of this manuscript and Fig. 9 
of Refill). From this point of view, the phase at 1/3 ML 
can be described by an ordered array of La-vacancies 
in the 1/2 ML structure. There is an effective repul- 
sion between La- vacancies due to the repulsion between 
La-atoms on neighboring A-sites. We describe the to- 
tal energy by an empirical model energy of the form 
E(cy) = E + Ef ■ cy + A ■ Cy, where cv is the con- 
centration of La vacancies, Eq is the energy of the struc- 
ture with all A-sites filled (1/2 ML), Ef is the formation 
energy of an isolated La-vacancy, and A describes the re- 
pulsion between vacancies. Coexistence between the two 
phases would result from a negative value of A. In that 
case, adding an additional ad-atom to a phase requires 
more energy than starting a new phase with the next 
higher coverage. Between 1/3 and 1/2 ML, however, A 
is positive as filling a portion of vacancies is favorable 
compared to creating patches of pure 1/2 ML coverage. 

We calculated the energy of an adsorption structure 
with three La atoms on neighboring A sites separated by 
one vacancy within one valley. La triplets in different 
valleys have been arranged, so that the distance between 
vacancies is maximized in order to minimize the repulsive 
energy. Based on the energies at 1/3 and 1/2 ML as well 
as at the intermediate coverage of 3/8 ML just described, 
we can determine the three parameters Eq, Ef and A to 
be 0.05, -0.56 and 0.26 eV, respectively. 

At a certain vacancy concentration of c v = 17% (i.e. 
a La-coverage of 0.42 ML) we find a phase boundary with 
the next stable phase at 2/3 ML at a chemical potential 
of 0.94 eV. According to our phase diagram, the pure 
surface reconstruction at 1/2 ML is never formed. The 
shaded region in Fig.^|corresponds to 1/2 ML structural 
template with variable vacancy concentration. 

As seen in the phase diagram shown in Fig.^Jbulk sili- 
cide formation becomes thermodynamically stable within 
the stability region of the 1/3 ML coverage. In a growth 
experiment we would expect the formation of bulk sili- 
cide grains to be delayed beyond a coverage of 1/3 ML. 
The nucleation of silicide grains may suffer from the large 
mismatch between bulk silicide phases and silicon. This 
is of particular importance during the initial stages of 
nucleation because the strained interface region occupies 
most of the volume of the grain. 

Thus it may be of interest to know the stability of 
silicide thin films on Si(001). We found one such silicide 
layer which is shown in Fig.^3 It consists of a (1 x 1) sil- 
icon surface in contact with two La layers that sandwich 
a layer of Si 4 ~ ions in between. While we have not per- 
formed a thorough search of other candidates, the energy 
of this silicide layer indicates that silicide formation will 
at the latest be initiated beyond a coverage of 2/3 ML. 
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FIG. 13: (Color online) Silicide overlayer at a coverage of 
2 ML. 




FIG. 14: Supercells used for the simulation of isolated La 
dimers shown in Fig. |3 



We can thus only pin down the onset of silicide formation 
within a coverage interval between 1/3 (thermodynami- 
cally) and 2/3 ML (including kinetic considerations). 



images introduced by the calculational supercell in or- 
der to emphasize that fact that this local arrangement 
corresponds to an isolated dimer. Fig. El sketches the 



a)<^° b) C) y~ d) 7^ 




FIG. 15: Supercells used for the calculation of La chains at 
1/6 ML as listed in Tab. EH 



supercells used. They were chosen in order to avoid frus- 
tration of Si dimers due to periodic images. 

Tab. ^0] summarizes the energetics of chains structures 
built from La dimers. The supercells used in the corre- 
sponding total energy calculations are sketched in Fig. 1151 



XI. CONCLUSIONS 



In this paper we have investigated the surface struc- 
tures of La adsorbed on Si(001) as a function of cover- 
age. We propose a theoretical phase diagram by relating 
the phase boundaries at zero temperature to chemical 
potentials, which can be converted into partial pressure 
and temperature in thermal equilibrium. 

Our findings elucidate the chemistry of third row ele- 
ments on Si(001) and the phases of La on Si(001), and are 
expected to provide critical information for the growth of 
a wide class of high-k oxides containing La. The phase di- 
agram may be used as a guide for the growth of La-based 
oxides on Si(001). 

APPENDIX A: SUPERCELLS AT DILUTE 
COVERAGES 

Fig. |21 shows the three possible La-dimer orientations 
on the Si(001) surface. We did not draw the periodic 
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